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Cosmological restrictions on conformally invariant SU{5) GUT models 
Abstract. Dirac's theory of constrained Hamiltonian systems is applied to the minimal 
conformally-invariant SU (5) grand-unified model studied at 1-loop level in a de Sitter 
universe. For this model, which represents a simple and interesting example of GUT 
theory and at the same time is a step towards theories with larger gauge group like SO{10), 
second-class constraints in the Euclidean-time regime exist. In particular, they enable one 
to prove that, to be consistent with the experimentally established electroweak standard 
model and with infiationary cosmology, the residual gauge-symmetry group of the early 
universe, during the whole de Sitter era, is bound to be SU{3) x SU{2) x U{1). Moreover, 
the numerical solution of the field equations subject to second-class constraints is obtained. 
This confirms the existence of a sufficiently long de Sitter phase of the early universe, in 
agreement with the initial assumptions. 



PACS numbers: 02.60.Cb, 04.20.Fy, 11.15.Ex, 98.80.Cq, 98.80.Dr 
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Cosmological restrictions on conformally invariant SU{h) GUT models 



1. Introduction 

In recent work by the first two authors (BucceUa et al 1992) the spontaneous-symmetry- 
breaking pattern of SU (5) gauge theory was studied in a de Sitter universe. The main result 
was the proof that the technique described in BucceUa et al 1980 to study spontaneous 
symmetry breaking of SU (n) for renormahzable polynomial potentials in flat spacetime can 
be generalized, for SU (5), to the curved background relevant to the inflationary cosmology 
(i.e. de Sitter). One thus obtained a better understanding of the result, previously found 
with a different numerical analysis (Allen 1985), predicting the slide of the inflationary 
universe into either the SU{^) x SU{2) x U{1) or SU{4) x U{1) extremum. 

The main tool used was the Wick-rotated path integral for Yang-Mills-Higgs theory at 
1-loop level about curved backgrounds (leading to a de Sitter model with S'^ topology), and 
the corresponding 1-loop effective potential F(r, ^) first derived in Allen 1985. Assuming 



of SU{5), and using the bare Lagrangian and the tree potential appearing in equations 
(2.1)-(2.2) of BucceUa et al 1992, the background-field method and the choice of 't Hooft's 
gauge-averaging term lead to (Allen 1985, BucceUa et al 1992) 




belongs to the adjoint representation 
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(1.1) 



167r2r4 



where Q has been defined in equation (2.6) of Buccella et al 1992, and we here denote by A 
the parameter defined in equation (2.7) of BucceUa et al 1992, to avoid confusion with the 
cosmological constant A. Note also that A denotes the special function defined in equation 



to the desired = 0) value. 

A naturally-occurring question is whether further restrictions on the SU{5) broken- 
symmetry phases can be derived within the framework of inflationary cosmology. This 
paper is devoted to the study of such a problem, and is thus organized as follows. Section 
2 performs the Hamiltonian analysis of the Riemannian (i.e. Euclidean-time) version of a 
de Sitter background coupled to the SU{5) model, following Dirac's theory of constrained 
systems. Section 3 studies the 1-loop effective potential in the SU{3) x SU{2) x U{1) 
and SU{4:) x U{1) broken-symmetry phases. Section 4 shows that to be compatible with 
the prescriptions of the electroweak standard model and with the inflationary scheme, the 
early universe can only reach the SU{3) x SU{2) x U{1) broken phase and that during 
all the de Sitter era, possible tunneling processes towards the SU{A) x U{1) invariant 
phase are energetically forbidden. Section 5 describes the numerical integration of the 
corresponding field equations for the four-sphere radius and the components of the Higgs 
field. Concluding remarks are presented in section 6. 



(A.l) of the appendix, r 



■\/3/A is the four-sphere radius, and the constant Vq is equal 
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2. Hamiltonian analysis 



The effective Lagrangian for the Riemannian version of an exact de Sitter background with 
topology coupled to the SU (5) model is given by (up to a multiplicative constant) 



4ttG 



(2.1) 



where the derivatives (pi are taken with respect to the Euclidean time r and the spatial 
gradient of the Higgs field has been assumed to be negligible. Note that, after having 
integrated over the gauge-bosons degrees of freedom, the effective potential only involves 
the Higgs field. According to the simplified argument usually presented in the literature, 
the field equations obtained by varying the action with respect to r and ^ are then (Allen 
1983) 

3 



8ttG 



E-=i /2 + V{r, *) 



SV{r, 4) 
6(pi 



\/i= 1,..,5 



(2.2a) 



(2.26) 



However, this approach does not take into account the full Hamiltonian treatment of the 
problem. In other words, from equation (2.1) one derives the primary constraint Pr ~ 0, 
where Pr is the momentum conjugate to the four-sphere radius, and ~ is the symbol of 
weak equality (i.e. an equality which only holds on the constraint surface). This primary 
constraint should be preserved using the technique described for example in Dirac 1964 
and Esposito 1992. The corresponding Hamiltonian analysis is as follows. 
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The canonical Hamiltonian He-, defined as the Legendre transform of the Lagrangian 
L in equation (2.1), takes the form 



5 2 

ffe = r"^ ^ ^ + - rV(r, *) (2.3) 

i=l 



where we have defined b = 3/(47rG). Thus, the effective Hamiltonian H defined on the 
whole phase space becomes 

H = Hc + X{r,^,Pr,p^) Pr . (2.4) 
The preservation of the primary constraint Vi = Pr yields 



XT/ 

Pr =[pr, ^ 2r-^ J^Pli - + 4rV(r, *) + r^— (2.5) 

i=l 



where | , | denote the Poisson brackets. One therefore finds in our model the secondary 



constraint 

V2 = 2r-^J2pl^-2br + 4:r^V{r,^)+r^^ . (2.6) 
Further constraints are not found, since the preservation of V2 leads to the condition 

= V2 ={V2, ^{V2, i^e} + X{V2, Fi} (2.7) 

which can be solved for A as 



V2.hA 

X = -^ f . (2.8) 

{V2,Vi\ 
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This solution can be obtained since the constraints Vi and V2 are second-class, and the 
Poisson brackets appearing in the formula for A are found to be 

8 A 5V 4^ <5V 

^=l %=1 



x2 T /" XT r 

1^2, Vi} = -10r-6^p2^ -26 + 12rV(r,*) + r4^ +8r3— . (2.10) 



Moreover, since V\ and V2 are second-class, they can be set strongly to zero using Dirac 
brackets (Dirac 1964, Esposito 1992), hereafter denoted by { , }*. The corresponding 
field equations are 

ff«|r,#}*RiA (2.11) 



{^^^HY^r-w{l-l^) (2.12) 



Pr^[pr,HY ^[pr,H^]-[pr.Vl^C^^[Vm.H^ ^ Q (2.13) 

where C/m is the matrix of Poisson brackets of second-class constraints. Note that, since 
V2 = when Dirac brackets are used, equations (2.12) and (2.14) can be written as 

~ r-'^Pv^ (2.15) 
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p ^ r^l^ . (2.16) 

For the purpose of numerical integration, the most convenient form of these equations is 



A (^4^^) ^^4^ (2.17) 



Such a system is here solved choosing the following initial conditions: 



(2.18) 



r(0) = ro = V3/87rGFo (2.19) 

Vm = (2.20) 

p^(0) = (2.21) 

P^,(0) = (2.22) 

where in equation (2.19) the value chosen for ro leads to a suitable cosmological constant for 
the inflationary era, and in equation (2.22) we have neglected for simplicity initial-kinetic- 
energy effects (cf end of section 5). Of course, the values should obey the constraint 
V2{t = 0) = 0, i.e. 



AT/ 

= -26ro + 4(ro) V(ro, *o) + (ro)^ ^ 



(2.23) 



Prom now on it is useful to use dimensionless units. For this purpose, we define a = 
^/2G/Stt = ^/2/(SttMp) and make the rescalings r^crr, r— >crT,<^i^ 0j/(7rcr-\/2), 
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and we also define 2TT'^a'^V{ar,(j)i/ {Tray/^)) = V{r,(f)i). The dimensionless field equations 
corresponding to equations (2.17)-(2.18) are then found to be 



+ ^^^^^ 

(-101-2 Eti - 2 + 12r2V + r<l0 + Sr^f^) 



Moreover, the initial conditions here chosen take the form 



r(0) = ro = (2.26) 



0,(0) = 0° (2.27) 

p^(0) = (2.28) 

p^,(0)=0 (2.29) 
where the 0° values obey, for a given ro value, the constraint 



= -2ro + 4(ro) V(ro, (/>°) + (ro)^ ^ 

or 



(2.30) 



Note that the choice Vq — fixes reasonably the critical temperature for the phase 
transitions to be of order Tc = Mx (Mx, approximately equal to 10"^^ Gev, is the typical 
order of magnitude of the unification mass in the minimal SU{5) model). 
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3. 1-loop effective potential 



We here study cj) in the forms invariant under the subgroups SU{2)) x SU{2) x U{1) and 
SU{A) X U{1) respectively, i.e. 



0321 = ^^^i^diag(2, 2, 2, -3, -3) 



30 



(3.1) 



11^ ||l/2 

.41 = ^^yUdiag(l, 1, 1, 1, -4) 



20 



(3.2) 



since it has been shown in AUen 1985 and BucceUa et al 1992 that these are the only 
subgroups relevant to the SUih) symmetry-breaking pattern. From now on, we denote by 
7(t) the norm of ^321 or of ^41, which is the only variable characterizing these broken- 
symmetry phases. 

The form of the dimensionless effective potential in the 321 and 41 directions (i.e. 
when (p — 0321 OT 4> — 04i) is obtained inserting equations (3.1)-(3.2) into equation (1.1), 
which yields 

V{r, 032i)|||^||i/2=^ = Vs2i{rn) 



45a f 1 
47rr2 1^+3 



1 - log 



37rM|> 



7 



+ 



25a" 
32^ 



^ , / 2M| 2 
1 - log! ^ ,5o r^ 



2r4"^ 



5a 2 2 
■^r 7 

DTT 



37rM|, 



+ Vo 



7 



(3.3) 
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V{r, 041 ) I II 1/2=^ = V4i{r,^) 



45a r 1 
47rr2 P + 3 



7 



+ 



75^2 1 



167r2 4 



X ^2 



A 

+ 3 ^7^ 



Set 
47r 



2 2 

r 7 



(3.4) 



where Vq = (8M_^)/(9M|,) and a = q^/Att. It is also useful to derive the r — > oo limit of 
these potentials (i.e. their flat-space limit) as 



~ , - , . \ ~ 25q: 

V'(</>32i)|||0||i/2=^(r ^ oo) = Fo + 32^7 



log 



4M|. 



■7 



(3.5) 



^(04i)|||<;||i/2^^(r ^ oo) = "l/o + 



■A7^ + "4 



167r2 ' 647r2 



7 



, /l5aM| o\ 1 

log 

^\ 8M 



2 ' 
X 



. (3.6) 



The potentials (3.5) and (3.6) evaluated at their minima, denoted for simplicity by 7^, 
turn out to be 

1 Mi 



^(032i)|||^l|i/2=^^ (r ^oo) = Vo- 



47r2 M| 



1 Mi 



^(<^4i)|||^l|i/2=^^(r ^ 00) = - ^]^^exp 



rA 



(3.7) 



(3.8) 



The experimental evidence for the SU{2>) x SU{2) x U{1) gauge symmetry at energy 
E greater than or of order 100 GeV requires for the scalar potential that the absolute 
minimum of V{(f)32i){f 00) should remain below the absolute minimum of V^(^4i)(r — > 
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oo) and below Vq. One thus finds the condition on the bare parameters A2 and A4 (Allen 
1985) of the scalar potential 



-A4-A,>-^log -) . (3.9) 



From now on, we restrict the choice of A2 and A4 to the region defined by the inequality 
(3.9). Under the assumptions described so far, the initial conditions (2.27) and (2.29) take 
the form 

7(0) = 70 (3.10) 

7(0) = . (3.11) 

The 7o value and the corresponding residual symmetry are obtained by solving separately 
the constraint equations (cf equation (2.23)) 



= -2ro + 4(ro)^F32i(ro,7o) + (ro) 



5r 



(3.12) 



ro,7o 



= -2ro + 4(ro)V4i(ro,7o) + (^o) 



5r 



(3.13) 



After doing this, one compares the V32i(ro,7o) and V4i(ro,7o) values, requiring that the 
correct initial condition should lead to the minimum value of the effective potential. Once 
the correct initial condition has been picked out in this way, the system (2.24)-(2.25) 
expressed in terms of 7 becomes 



— (r^) f« r^— (3.14) 
ctr ^ ^ ()7 
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(3.15) 



jointly with the constraint (cf equation (2.6)) 



5r 



~ AT/ 

V2{t) = 2r^^^ - 2r + 4rV(r, -f) + r^— (3.16) 

or 



which should vanish Vr. 



4. Absolute minimum 



By virtue of equations (A.7)-(A.19) of the appendix, the asymptotic form of the constraints 
(3.12)-(3.13) is 



25a 
32^ 



70^ 



, fhaMl 2 



X 



+ 



(^0-1 







25^2 
32^ 



7o 



3, /15aM2 2 



3 6~ 

4 + 5^ 



+ O(ro-2)=0 



(4.1) 



(4.2) 



The numerical solution of equation (4.1) yields 70 = 1.42 • 10~^, where the NAG-library 
routine C05ADF has been used in double-precision. Such a value of 70 is compatible with 
the asymptotic formulae appearing in the appendix. Interestingly, if the inequality (3.9) 
holds, one finds Vr 



^4i(t)-F32i(t) 



25^2 
32^ 



7 



3, /3\ 1 6- 
2^"<2)-4+5^ 



> 



25^2 
32^ 



7 



3 /3x 1 
4^"<2j-4 



>0 . (4.3) 
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This result ensures that at r = 0, if tq is taken to be of order 10®, so that the asymptotic 
formulae of the appendix can be applied, the only possible residual symmetry is SU (3) x 
SU{2) X U{1). Remarkably, because equation (4.1) is, at least up to first order in the 
ro-expansion, independent of the bare parameters (A, Q) of the scalar potential and only 
dependent on Mx, the value found for 70 is basically model-independent. If the alternative 
symmetry SU{4:) x U{1) were chosen, the inequality (4.3) would lead to a tunneling of the 
Higgs field towards the energetically more favourable phase. Moreover, our result (4.3) 
ensures that, once the initial SU{3) x SU{2) x ^7(1) symmetry is chosen, the Higgs field 
remains in this broken-symmetry phase during the whole de Sitter phase of the early 
universe. 

5. Numerical analysis 

The numerical integration of the system (3. 14)- (3. 15) can be performed after reduction to 
first-order form, and using the second-class constraint (3.16) and equations (2.15)-(2.16). 
Thus, defining 



7(t) = t]{t) 



(5.1) 



this leads to 




SV_ 
5^ 




(5.2) 
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r{T)^-, ^ '-2^ ^ . (5.3) 

M2 - 32rW - 13r3^ - r^^) 

By taking ro = 10^, which is a typical value for GUT models in a de Sitter universe, 
the system (5.1)-(5.3) has been solved within the framework of the Runge-Kutta-Merson 
method (NAG- library routine D02BAF). The results of our numerical analysis are shown 
in figures 1 — 3, which present the Euclidean-time evolution of 7, 7 and r respectively. 
These plots clearly show the existence of an almost exact de Sitter phase (i.e. when the 
four-sphere radius remains approximately constant) with typical time r, in dimensionless 
units, of order 3 • 10^, whereas the duration of the exponentially-expanding phase may 
be taken to be of order 7 • 10^. Interestingly, this result is practically independent of the 
bare parameters in the SU (5) potential we have chosen if the inequality (3.9) is satisfied. 
Further restrictions are then given by particle physics, i.e. proton-lifetime experiments and 
renormalization-group equations for the coupling constants and (smOw)^ (Buccella et al 
1989, Becker-Szendy et al 1990, Amaldi et al 1992). 

If T > 7 • 10^, figure 3 shows a rapid variation of the four-sphere radius, so that the 
early universe is no longer well described by a de Sitter or exponentially-expanding model. 
Note also that the numerical analysis here presented rules out the occurrence of tunneling 
effects between broken-symmetry phases. Moreover, it should be emphasized that, setting 
to zero the initial kinetic energy (cf equation (3.11)) one obtains the most favourable initial 
conditions for a long de Sitter phase of the early universe, whereas values of 7(0) 7^ may 
be shown to lead to a much more rapid variation of the four-sphere radius r. 
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6. Concluding remarks 

This paper has studied the SU (5) symmetry-breaking pattern in a de Sitter universe from 
a Hamiltonian point of view. As a result of this analysis, one finds that the model is 
characterized by two second-class constraints. Interestingly, the secondary second-class 
constraint has been used to prove that the early universe can only reach the SU{2>) x 
SU (2) X broken-symmetry phase if we require the correct low-energy limit for the GUT 
theory. Furthermore, under this obvious requirement, the SU{2>) x SU{2) x U{1) invariant 
direction turns out to be energetically more favourable than SU{A) x U{1). Thus, during 
the whole inflationary era, possible tunneling effects between the two broken-symmetry 
phases are forbidden. 

This conclusion supersedes earlier work on the SU (5) symmetry-breaking pattern in 
de Sitter cosmologies appearing in Allen 1985 and Buccella et al 1992. It also provides 
a relevant example of cosmological restrictions on GUT models (cf Collins and Langbein 
1992). 

Moreover, the resulting field equations have been solved numerically. Standard meth- 
ods for the numerical integration of first-order systems of ordinary differential equations 
show that the early universe starts out in a de Sitter state. The exact dc Sitter state is 
then replaced by a more general exponentially-expanding universe corresponding to the 
slow-rolling-over phase, as one would expect. In dimensionless units, the total duration of 
these two phases has been found to be of order 7 • 10^. 
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From a field-theoretical point of view, the 1-loop effective potential of equation (1.1), 
originally derived in Allen 1985, and also used in Buccella et al 1992 and in this paper, 
might be subject to criticisms, since topological methods have been used to prove the 
general non-existence of a gauge choice for a theory with gauge group SU (N), N >3, over 
the manifolds and S"^ (Jungman 1992 and references therein). However, this difficulty 
holds for a classical non-Abelian gauge-field theory. By contrast, we have studied the 
quantization of a Yang-Mills-Higgs theory where gauge-averaging terms are included in 
the Wick-rotated path integral using the Faddeev-Popov technique and a specific choice 
first proposed by 't Hooft (Allen 1985, Buccella et al 1992). Although the choice of a 
gauge-averaging term is clearly suggested by what would be set to zero at the classical 
level, there is by now some evidence that the resulting quantum theory is not equivalent 
to the model where only physical degrees of freedom are quantized after fixing the gauge 
(Esposito 1992 and references therein). It therefore appears that the approach used in Allen 
1985 to evaluate the 1-loop effective potential for SU (5) theory in de Sitter cosmologies 
remains a useful tool to improve our understanding of the symmetry-breaking pattern of 
non-Abelian gauge theories in curved backgrounds. 

Appendix 



The function A{z) appearing in equation (1.1), and its derivatives, are given by (Allen 
1985) 
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'1 



- y{y -^){y -^)'^{y) dy 
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A'{z) = --{z + 2) 



z 1 
+ 2 + 3 



(^.2) 



A"{z) 



1 1 

2 ~ 2 



, 1 + 
4 /Ti 



(-4.3) 



At large the following asymptotic expansions hold (Allen 1985): 



( z^ 19\ 3 

- ( ^ + ^ + ^ ) log (^) + g + ^ + const. ^0{z-^) (A.4) 



^'(^)~-l(^ + 2)log(^) + | + 0(^-^) 



M\z)^-\\o^{z)^0[z-^) . 



(A.6) 



In section 4, we rely on the following exact formulae for the derivatives of the dimensionless 
effective potential in the SU{3) x SU{2) x U{1) and SU{4:) x U{1) broken-symmetry phases: 
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In light of equations (A.4)-(A.6), the asymptotic expansions of the dimensionless effective 
potential and of equations (A.7)-(A.12) are given by 
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Figure captions 



Figure 1. The Euclidean-time evolution of the norm 7 of the Higgs field is here plotted, 
after solving numerically the system (5.1)-(5.3). 



Figure 2. The evolution of 7 is here shown. 



Figure 3. The numerical solution for the Euclidean-time evolution of the four-sphere 



radius r is here presented. If r e 



0,3 • 10^ 



, our solution of the system (5.1)-(5.3) 



shows that r is approximately constant, as one would expect during the de Sitter 
phase of the early universe. 
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